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Abstract. We study the influence of scale-free correlated disorder on the metal-insulator transition in 
the Anderson model of localization. We use standard transfer matrix calculations and perform finite-size 
scaling of the largest inverse Lyapunov exponent to obtain the localization length for respective 3D tight- 
binding systems. The density of states is obtained from the full spectrum of eigenenergies of the Anderson 
Hamiltonian. We discuss the phase diagram of the metal-insulator transition and the influence of the 
correlated disorder on the critical exponents. 
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1 Introduction 

The possibility of having phase transitions in disordered 
systems, which contain randomness as a central ingredi- 
ent, has attracted a lot of interest over past decades. Pro- 
totypical examples of classical and quantum disordered 
systems are the percolation problem [l] and the Anderson 
model of localization [2] , respectively. Over the years, the 
respective transitions - percolation and Anderson metal- 
insulator transition (MIT) - have been intensively stud- 
ied [3j|7] . Typically, the random numbers, which represent 
the disorder, are taken to be uncorrelated. However, in 
realistic systems, where, for example, the disorder is in- 
duced by a complex environment surrounding the system 
sites, one expects to find correlations between the ran- 
dom numbers. Spatial correlations can be characterized 
according to their behavior on different length scales. For 
example, they might become irrelevant if the length scales 
associated with the phase transition are larger than a char- 
acteristic correlation length. On the other hand, there is 
also scale-free disorder, which is found in many physical 
systems j5j[Sj[0] . Here, the correlations are taking effect on 
all length scales. These long-range correlations are char- 
acterized by a power-law behavior, C(r — r') cx |r — r'| _Q . 
Here, C(r — r') denotes the correlation function and a is 
the correlation exponent. 

For the (classical) percolation problem it was found 
that the presence of scale-free disorder has a profound in- 
fluence on the percolation transition. For this situation, 
the extended Harris criterion [10}fl3] predicts a crossover 
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of the critical exponent v from its value i>q for uncorre- 
lated random numbers to 2/ a provided the decay of the 
correlations is sufficiently weak, a < 2/vq. 

In the present paper we address the question of the 
role of power-law correlations for the Anderson transi- 
tion in disordered electronic systems. Originally, in his 
seminal paper Anderson showed that extended electronic 
states can become spatially localized due to the presence 
of uncorrelated disorder [21. As a consequence, the sys- 
tem undergoes a phase transition from a conducting phase 
(for extended states) to an insulating phase (for localized 
states), which can be characterized by a critical exponent. 
In a previous study [14] it was found, that the critical 
exponent is independent of the correlation exponent for 
a transition at fixed energy in the center of the band, 
while for a transition at fixed disorder strength the crit- 
ical exponent obeys the extended Harris criterion. How- 
ever, the calculations have been performed using a mod- 
ified transfer-matrix method (TMM), which consists of 
forward and backward TMM calculations for a quasi-one- 
dimensional (quasi-lD) block of length Lq ~ 10 3 . This 
artificial periodicity might have an additional influence 
on the transition. To avoid this issue, we use the stan- 
dard TMM [7] for calculating the localization length of 
quasi-lD systems with a length L = 4 • 10 5 . Subsequently 
performing a finite-size scaling (FSS) analysis provides us 
with estimates of the critical points [15] for different cor- 
relation exponents. The critical points are summarized in 
a phase diagram showing the influence of the correlations 
on the phase boundary. This is the central result of the 
present paper. Additionally, we calculate the density of 
states (DOS) of 3D systems in the presence of scale-free 
disorder. 
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The paper is organized as follows. In the next sec- 
tion, we introduce the Anderson model of localization and 
briefly summarize the main properties of the associated 
phase transition. Moreover, we provide an overview of the 
numerical methods we use to calculate the properties of 
the transition. In Sec. [3] we present our results for transi- 
tions at fixed energy and fixed disorder. Then the phase- 
diagram of the Anderson MIT in the presence of scale-free 
disorder is discussed and the relation to the DOS is in- 
vestigated. Finally, in the last section we summarize and 
discuss our results. 



2 Model and Numerical Methods 

2.1 Anderson Model of Localization with Long-Range 
Correlated Disorder 

The Anderson model [2]7] is widely used to investigate the 
phenomenon of localization in disordered materials. It is 
based upon a tight-binding Hamiltonian in site represen- 
tation 

H = 2>|i><i|-x> u |i><j|, a) 

i ij 

where |i) is a localized state at lattice site i. The ma- 
trix elements iy denote hopping integrals between states 
at sites i and j. Typically, hopping is restricted to near- 
est neighbors. The on-site potentials Si are random num- 
bers, chosen according to some probability distribution 
P{e) characterized by the mean (ej) and the correlation 
function C(|-£|) cx X]i( £ i £ i+^)- However, usually the site 
energies are taken to be statistically independent. For ex- 
ample, convenient choices of P(e) are a box distribution 
of width W or a Gaussian white noise distribution, both 



t2 Wo- 



Other distributions 



with (si) = and C(\£\) 
have also been considered [7j 16 17 

For uncorrelated potentials the resulting situation may 
be summarized as follows Vf\: for strong enough disorder, 
W > W c (0), all states are exponentially localized to a 
region of finite size. The extent of this region is charac- 
terized by the so-called localization length A. The value of 
the critical disorder strength W c depends on the distribu- 
tion P(e) and the dimension d of the system. The value 
of W c additionally depends on the Fermi energy E and 
the curve W C (E) separates localized states, W > W C (E), 
from extended states, W < W C (E), in the phase diagram. 
If instead of E the disorder strength is fixed, there will be 
a critical energy E C (W) and states with \E\ < E c are ex- 
tended and those with \E\ > E c localized. The transition 
from extended to localized wave-functions at the critical 
point is called disorder driven or Anderson MIT. In the 
vicinity of the critical point the localization length behaves 
as 

A(r) cx |r c - t\~ v , (2) 

where r is either E or W. The critical exponent v charac- 
terizes the phase transition and is expected to be univer- 
sal. 



In the present work we are interested in the influence 
of long-range correlated disorder potentials on the Ander- 
son MIT. In particular, we study the dependence of the 
critical points and the critical exponents on the strength 
of the correlations. To this end we use random potentials 
generated from a Gaussian probability distribution and 
with a correlation function of the form 



C{£) ee ( £i e i+e ) cx \£\- 



(3) 



where a is the correlation exponent which determines the 
strength of the correlations. In contrast to short-range 
correlations the power-law behavior in Eq. Efy does not 
introduce a characteristic length scale and therefore the 
disorder is said to be scale-free. 

In general, it is extremely complicated to obtain an- 
alytical results of transport properties for the Anderson 
model of localization. For example, only in the case of 
d = 1 rigorous proofs of strong localization for all en- 
ergies and disorder strengths have been given 18]. More- 
over, the explicit energy and disorder strength dependence 
of the localization length for weak disorder has been de- 
rived [19 20 . There are also some results for ID systems 
with long-range correlated disorder. For energies close to 
the band center a weak disorder expansion for A has been 
derived in Ref. [21], which shows the dependence of the 
localization length on the correlations via a Fourier trans- 
form of the correlation function. Here, the localization 
length is proportional to W~ 2 . On the other hand at the 
unperturbed band edge (\E\ = 2) it was found 
that X(E = 2, W) oc W~ v with y = 2/3 for a = 
y = 2/(4 -a) for a < 1. 



2.2 Numerical Methods 

In order to investigate the influence of scale-free correla- 
tions on the Anderson MIT, we first generate the corre- 
lated on-site potential for systems of size M x M x L us- 
ing a modified Fourier filtering method (FFM) [24] with 
one additional step. Namely, after performing the usual 
FFM we shift and scale the obtained sequence of corre- 
lated random numbers such that the mean vanishes and 
the variance is W 2 /12. 

As mentioned in the introduction, the localization length 
A is calculated using a standard TMM m. Thereby we 
use a new seed for each parameter combination (E, W, 
a, M). Lastly, the critical exponent, mobility edge and 
critical disorder are obtained from the FSS analysis [15] 
based on a higher-order expansion of Eq. ([2| at the tran- 
sition (r = t c ). This procedure is outlined in appendix [A] 
The error resulting from the associated fitting procedure 
should not be seen as an upper (or lower) bound for the 
critical parameters, but rather as a qualitative measure for 
the stability of the fitting procedure. A reliable estimate 
of the numerical uncertainty requires a more sophisticated 
error analysis [25] . 

Another issue connected with the FSS method are the 
corrections to scaling. Generally, for small systems one 
would always expect to find finite-size corrections and con- 
sequently one should include corrections to scaling in the 
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analysis. However, this also increases the number of pa- 
rameters to be fitted tremendously and makes it some- 
times complicated to find a reasonable fit. One possibility 
to partly circumvent this problem is ignoring the small 
system sizes (in our analysis below this means M = 5,7) 
and doing the FSS without corrections. 

The DOS is obtained from the full spectrum of eigenen- 
ergies of the 3D Anderson Hamiltonian for systems of 
size M 3 . The eigenenergies are calculated using standard 
matrix diagonalization methods 26 . Since the accessible 
maximum system size is restrictea by the numerical re- 
sources, the results are averaged over a large number of 
disorder realizations to decrease statistical fluctuations. 
Furthermore the symmetry of the DOS with respect to 
E = is utilized. 



3 Numerical Results and Discussion 
3.1 Numerical Calculations 

In order to study the localization length in the presence 
of correlated disorder we focus on quasi- ID systems with 
L = 400000 and M = 5,7,9,11,13 and 15. The error 
of the localization length is determined from the vari- 
ance of the change of the Lyapunov exponent during the 

The accuracy of the localization 



27 28 



TMM iterations ^ 
length is thereforehmited by the finite length of the con- 
sidered systems. To give an impression of the quality of 
the TMM results and the FSS fitting procedure Fig. [I] 
shows the reduced localization length A = X/M versus 
disorder strength W for E = 6.0 and a — 2.5. In Fig. [lji 
the raw data for A obtained from the TMM calculation 
are shown for various system sizes. Performing the FSS 
procedure taking corrections to scaling into account one 
obtains the curves shown in Fig. [T^,. Considering no cor- 
rections to scaling we obtain the fits shown in Fig. [T]d. In 
this case, ignoring the small systems leads to almost the 
same critical values. In the following we will concentrate 
on results obtained without taking corrections to scaling 
into account. 

The DOS is computed for disorder strengths W — 1.5, 
. . . , 30 for the uncorrelated and a long-ranged correlated 
potential (a — 0.9), respectively. The size of the systems 
is M 3 = 22 3 . Results are averaged over at least 1000 disor- 
dered samples. For the ordered system, W — 0, the DOS 
is calculated by the diagonalization of a single system with 
M 3 = 30 3 . 



3.2 Transition at Fixed Energy 

First we focus on the Anderson MIT for fixed energy 
(E = const.). The values of the respective critical param- 
eters are shown in Tab. [JJ One sees that for uncorrelated 
disorder (a = oo) the obtained critical exponents, v, are 
consistent with the high-precision value vq — 1.58 ± 0.03 
of Ref. 15 . For E = the critical disorder strength 
agrees very well with the value found previously [TBI , W c = 
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Fig. 1. Reduced localization length A vs disorder strength 
W for E = 6.0 and a = 1.5. Solid lines show the FSS fit 
to numerical data (a) for all system sizes M taking correc- 
tions to scaling into account and (b) for M > 7 without 
corrections to scaling. The resulting critical exponent v 
and disorder strength W c are shown together with the re- 
spective expansion orders of the FSS (cf. appendix |A|) . 



In the presence of long-range correlations the critical 
exponent remains close to the value for uncorrelated dis- 
order potentials at least for energies well inside the band 
of the system without disorder. For energies close to the 
band edge (\E\ — 6) of such systems the critical expo- 
nent is smaller than vq in case of a long-range correlated 
potential. However, at the same time the estimated error 
of the exponents becomes larger and is more sensitive to 
changing the fitting parameters. The error for the critical 
disorder is relatively small independently of the value of a 
and one finds that the value is more stable than v against 
changing the fitting parameters. 

The value of W c is monotonically increasing for de- 
creasing a. In other words the MIT sets in for a larger 
disorder strength compared to the uncorrelated case. This 
supports the intuitive expectation of an effective smooth- 
ening of the random potential due to the correlations. For 
ID systems and weak disorder it has been shown that the 
effective disorder is given in terms of the Fourier transform 
of the correlation function [21 . In the center of the band 
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this leads to an increasing localization length for smaller 
correlation exponents (29). 



3.3 Transition at Fixed Disorder Strength 

Next, we consider the case of transitions at fixed disorder 
strength (W = const.). The respective results for criti- 
cal energies and exponents are summarized in Tab. [2] For 
comparison Tab. [3] contains some FSS results obtained 
with corrections to scaling taken into account. 

For uncorrelated random potentials previous studies 
showed that the critical exponent of the transition at fixed 
disorder is, within error bars, identical to the exponent ob- 
tained for the transition at fixed energy. Here, we also find 
a good agreement with the respective exponent shown in 
Tab. [I] and with the high- precision value v$ [15]. Further, 
the critical disorder strengths are in accordance with the 
results of Ref. 



30 



In contrast to the transition at fixed energy discussed 
earlier, in the presence of long-range correlations the crit- 
ical exponents are larger than the respective exponent ob- 
tained for uncorrelated potentials for all disorder strengths 
except W — 6.0. Moreover, the mobility edge is systemat- 
ically shifted towards higher energies. 



3.4 Phase Diagram and DOS 

Combining Tabs. [T] and [2] we obtain a complete phase di- 
agram of the Anderson model in presence of long-range 
correlated disorder, which is shown in Fig. [2] The phase 
diagram reflects the general features we have discussed for 
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Table 1. Critical disorder W c and exponent v obtained 
from FSS analysis without taking corrections to scaling 
into account. The errors indicate the confidence interval of 
the fit. The symbols f and % denote parameters coinciding 
in Tables [Hand [2l 



the two transitions. In the presence of long-range correla- 
tions the metallic phase space grows, pushing the mobility 
edge to larger disorder strengths and higher energies. 

Figure [3] illustrates the influence of long-range corre- 
lated disorder on the DOS. The contours in Fig.[3]show the 
characteristic broadening of the DOS for increasing disor- 
der strength W [30] , The difference between correlated 
and uncorrelated disorder is much less pronounced than 
in the phase diagram. Around the band center (E < 2) the 
DOS is increased by correlations. Toward the band edges 
the DOS is slightly smaller. From Figs. [2] and [3] we can 
conclude that the mobility edge is always clearly inside 
the band and comes close to the band edge only for small 
disorder strengths W. 

One might also ask how the two transitions behave at 
the same point in the phase diagram. A previous study of 
the Anderson model with uncorrelated disorder suggested 
that close to the band edge and for fixed disorder strength 
the critical exponent may be different from v§ [31] . How- 
ever, strong finite-size effects in this region did not allow 
a conclusive answer. A more recent study indicates that 
for both transitions at different points in the phase dia- 
gram the same exponent is obtained [32] . From Fig. [2] we 
see that there are two pairs (Wc, £ , c jwhich denote the 
same phase-diagram point, respectively. The critical pa- 
rameters for these points are marked in Tabs. [T] and [2] 
by f and %. In Fig. [2] these points are conspicuous, be- 
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Table 2. Critical energy E c and exponent v obtained from 
FSS analysis without taking corrections to scaling into 
account. The errors indicate the confidence interval of the 
fit. 
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Table 3. Critical energy E c and exponent v from FSS 
analysis with taking corrections to scaling into account. 
The errors indicate the confidence interval of the fit. 
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Fig. 2. Phase diagram of the Anderson model with 
long-range correlated disorder. Open and filled symbols 
indicate transitions at fixed energy and fixed disorder 
strength, respectively. Lines are a guide to the eye only. 



transitions at fixed disorder strength are larger than vq. 
A similar behavior has been reported in Ref. [14] , where 
in addition an agreement of the critical exponents of the 
transition at fixed disorder strength with the extended 
Harris criterion was found. Due to the limited accuracy of 
our numerical data and the sensitivity of the critical expo- 
nent, we cannot give a quantitative comparison with the 
Harris criterion. Nevertheless, it is interesting to notice, 
that in our case the critical exponents strongly depend on 
the value of W and only weakly on the value of a. More- 
over, there is a possible connection of the behavior of v 
and the slope of the curve describing the phase boundary, 
W C (E) and E C (W), respectively. For example, for fixed 
disorder strength the critical exponents increase with in- 
creasing magnitude of dE c (W) / dW . In other words, only 
if the chosen path through the phase diagram is perpen- 
dicular to the phase boundary, we obtain an unchanged 
critical exponent from the FSS analysis compared to the 
uncorrelated value Uq. Otherwise the estimated exponent 
is different from vq. A detailed investigation of this behav- 
ior would certainly be very interesting and might help to 
elucidate the role of long-range correlations for the Ander- 
son MIT. 



cause the (red) open squares and filled circles coincide. In 
both cases the transitions at fixed energy yield a critical 
exponent in agreement with z/q, while the exponents of 
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Fig. 3. The contour diagram shows the disorder depen- 
dence of the DOS for the uncorrelated (black) and a cor- 
related potential with a — 0.9 (red dashed). Lines are 
drawn at the given percentage of the maximum value of 
the DOS located at W = and E = 0. The band edge 
(BE) is estimated by the largest eigenenergy observed. 



4 Summary and Conclusions 



In summary, we have studied the role of scale- free disorder 
in the Anderson MIT. The correlations are characterized 
by a power-law with a correlation exponent. The charac- 
teristics of the Anderson transition have been obtained 
from the numerically calculated behavior of the localiza- 
tion length in quasi-lD systems. We employed a standard 
TMM computation and estimated the critical exponents 
and critical points using a FSS analysis for different corre- 
lation exponents. Further, we obtained the phase diagram 
for the Anderson MIT in presence of scale-free disorder. 

We observe a shift of the phase boundary towards 
higher energies and stronger disorder, respectively. The 
latter may be understood as a result of an effective smooth- 
ening of the disorder potential in presence of correlations. 
A similar behavior has been observed for ID systems, 
where the localization length increases for smaller corre- 
lation exponents. 

Regarding the critical exponents we cannot draw quan- 
titative conclusions due to the high sensitivity of the fitted 
results. However, qualitatively we see strong indications 
that the critical exponents behave differently for transi- 
tions at fixed energy and fixed disorder strength as it was 
reported before 14 . For fixed energies \E\ > 6 the critical 



exponent remains consistent with the value for uncorre- 
lated disorder, while for fixed disorder strengths W the 
exponent increases for increasing W . Further investiga- 
tions in this direction would certainly be helpful to get 
a better understanding of the role of correlations for the 
Anderson MIT. 
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A Finite-Size Scaling 

A problem one is always faced with when using numerical 
methods to investigate phase transitions, is the fact that 
for finite systems there can be no singularities induced 
by a transition and the divergences are always rounded 
off ~ 



33 



However, the phase transition can still be stud- 
ied using FSS. Specifically, near the MIT one expects the 
following one-parameter scaling law for the reduced local- 



ization length 33 



M 



A{M,T,b)=F[~,x{T)b 1 ! v ,<t>{T)b-* 



(A.l) 



where b is the scale factor, % is a relevant scaling variable, 
<j> is an irrelevant scaling variable, v > is the critical 
exponent and y > is the irrelevant scaling exponent. 
The irrelevant scaling variable allows us to take account 
of corrections to scaling due to the finite size of the sam- 
ple. Here, the parameter r measures the distance from the 
mobility edge E c , r = \E — E c \/E c , or the distance from 
the critical disorder strength W c , r = \W — W c \/W c . The 
choice b — M leads to the standard scaling form 

A(M, t) = F{M 1 ' v X {t), M-^(t)) (A.2) 

with F being trivially related to J- '. For r close to zero we 
expand F into a Taylor series up to order n\ and obtain a 
series of functions F n [15 



A(M,t) = ]T rM^FnixM 1 ^) 



(A.3) 



n=0 



Each function F n is then expanded up to order tt-r . Ad- 
ditionally, x an d 4> are expanded in terms of the small 
parameter r up to order tor, and mi, respectively. This 
procedure gives 



(A.4) 



ri=0 



From Eqs. (A.2 1 and (A.3 1 one can see that a finite system 



size leads to a systematic shift of A with M, where the 
direction of the shift depends on the boundary conditions 



33]. Consequently, the curves A(M,t) do not intersect at 
the critical point r = for different system sizes. The term 
Fq on the other hand shows the expected behavior. Using 
a least squares fit of the numerical data to Eqs. (A.3) and 
(A.4) allows us to extract the critical parameters v, E c 
For the actual orders of the expansions 



and W c 15 17 



as given in the legends of Figs. [T^, and [T]d, we have to 
determine, respectively, 8 and 4 independent combinations 
of the expansion coefficients. 
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